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New Evaluation of Proton Structure 
Corrections to Hydrogen Hyperfine 
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Abstract: We consider the proton structure corrections to hydrogen ground-state hyperfine 
structure, focusing on a state-of-the-art evaluation of the inelastic nucleon corrections — the 
polarizability corrections — using analytic fits to the most recent data. We find a value for 
the fractional correction A po i of 1.3 ± 0.3 ppm. This is 1-2 ppm smaller than the value of 
Ap i one would deduce using hyperfine splitting data and elastic proton structure corrections 
obtained from modern form factor fits. In addition, we discuss the derivations of the relevant 
formulas, paying attention to lepton mass effects and to questions surrounding the use of 
unsubtracted dispersion relations. 

PACS Nos.: 31.30.Gs, 32.10.Fn, 14.20.Dh, and 13.40.Gp 
Resume : 

French version of abstract (supplied by CJP) 
[Traduit par la redaction] 

1. Introduction 

Hyperfine splitting (hfs) in the hydrogen ground state is measured to 13 significant figures in fre- 
quency units [1], 

£hfs(e~p) = 1 420.405 751 766 7(9) MHz . (1) 

This accuracy is especially remarkable to theorists, who are hopeful of obtaining calculations accurate 
to a part per million (ppm) or so. We are close to reaching this goal, but some improvement is still 
needed and there currently seems to be a few ppm discrepancy between the best calculations and the 
data. 

The main uncertainty in calculating the hfs in hydrogen comes from the hadronic, or proton struc- 
ture, corrections. Calculation in this case is not ab initio calculation, which is currently impossible, but 
rather calculations that relate proton structure contributions in the bound state energy to information 
about proton structure obtained from elastic and inelastic electron-proton scattering measurements. 
Historically, the elastic and inelastic contributions, the latter also called polarizability corrections have 
often been treated separately, with the elastic corrections further divided into a nonrelativistic Zemach 
term and relativistic recoil corrections. From a modern viewpoint, the elastic and inelastic corrections 
should be treated as a unit since the sum lacks certain ambiguities that exist in the individual terms. 
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This talk will focus on the polarizability contributions, but noting the last remark, discussion of the 
Zemach and recoil corrections will not be omitted. 

2. Hyperfine splitting calculations 

The calculated hyperfine splitting in hydrogen is [1, 2, 3] 

E h[s (e-p) = (1 + Aqed + A^ cak + Astr) K , (2) 
where the Fermi energy is 

37r 6m e m p 

with m r = m e m p /(m p + m e ) being the reduced mass (and there are hadronic and muonic vacuum 
polarization terms [2] which are included as higher order corrections to the Zemach term below). The 
QED terms are accurately calculated and well known. They will not be discussed, except to mention 
that could be obtained without calculation. The QED corrections are the same for muonium, so it is 
possible to obtain them to an accuracy more than adequate for the present purpose using muonium hfs 
data and a judicious subtraction [4, 5]. The weak interaction corrections [6] also will not be discussed, 
and are in any case quite small. We will discuss the proton structure dependent corrections, 

Astr = A e; + A inel = A Z + A P R + Ap i , (4) 

where the terms on the right-hand-side are the Zemach, recoil, and polarizability corrections. 

The proton structure corrections come from two-photon exchange, as diagramed in Fig. 1. Gener- 
ically, the diagram consists of Compton scattering of off-shell photons from an electron knit together 
with similar Compton scattering from the proton. (We are neglecting the characteristic momentum 
of the bound electron. This allows a noticeably simpler two-photon calculation than for a scatter- 
ing process [7]. One can show that keeping the characteristic momentum would give corrections of 
0(am e /m p ) smaller than terms that are kept [8].) 




Fig. 1. Generic two-photon exchange diagrams, giving proton-structure corrections to hyperfine splitting. 



If the hadronic state represented as a blob in Fig. 1 is itself just a proton, the two-photon exchange 
specializes to Fig. 2. The photon electron vertex is known, and we can calculate the diagram using [9, 
10] 

T ll = llJl F 1 {q 2 ) + ^-a^q v F 2 {q 2 ) (5) 

for the photon-proton vertex with incoming photon momentum q. The functions F\ and F 2 are the 
Dirac and Pauli form factors of the proton. They are measured in elastic electron-proton scattering, and 
we shall refer to the contribution to the hyperfine energy from Fig. 2 as the elastic contribution. 

One should know that using Eq. (5) is an approximation, because this form of the vertex is sufficient 
and correct only when both the proton entering and leaving the vertex are on mass-shell. However, if 
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Fig. 2. Two-photon exchange diagrams for the "elastic" proton-structure corrections to hyperfine splitting. 



one is able to do the box diagram calculation using an unsubtracted dispersion relation, the on-shell 
parts of the elastic calculation are all that is needed. The story fits nicely into the discussion surrounding 
the derivation of the inelastic contributions, and we reserve its telling until we reach that discussion. 
For now, we will just note that the final sum of the inelastic and elastic proton structure corrections 
includes terms precisely as derived from the above form of the photon-proton vertex. 
The elastic contributions are separated as 



F 2l 

Ep 



(6) 



where the separation is into non-relativistic and relativistic terms — "Zemach" and "recoil." Non-relativistic 
means the limit m p — > oo with m e held fixed and with proton size held fixed. (Proton size information 
is embedded in the form factors F\ and F 2 .) 

The Zemach correction was worked out by Zemach in 1956 [11]; in modern form it is 



A z = 



8am r 



dQ 



G E {-Q Z ) 



2 ,G M (-Q 2 ) 



1 



1 



-2am r rz 



(7) 



the last equality defining the Zemach radius rz and we have used Q 2 = —q 2 . The charge and magnetic 
form factors are linear combinations of f\ and F 2 , 



G 



M 



Ge 



F 1 +F 2 
n 

Ft 



Am 2 



F 2 . 



(8) 



The relativistic elastic corrections A^,, also known as recoil corrections, won't be displayed in this 
talk because they are somewhat long (although not awful; see [9, 10]). An important point about them 
is that they depend on the form factors and hence upon the proton structure. However, evaluating them 
with different analytic form factors based on fits to the scattering data reveals them to be decently well 
determined by present standards. 

When the blob in Fig. 1 is not a lone proton, we obtain inelastic contributions or polarizability 
contributions [12, 13, 14, 15, 16, 17]. The inelastic contributions are not calculable ab initio. Instead, 
one relates them to the amplitude for forward Compton scattering of off-shell photons off protons, 
given in terms of the matrix element 



T^(q,p, S) = ^— J d 4 C (pS\ j ^(OMO) \pS) 



(9) 



where is the electromagnetic current and the states are proton states of momentum p and spin 4- 
vector S. The spin dependence is in the antisymmetric part 



{H^u, a 2 ) + H 2 (u, a 2 )) - H 2 (u, q 2 )- 



S-qp" 



p-q 



(10) 
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There are two structure functions Hi and H 2 which depend on q 2 and on the photon energy v, defined 
in the lab frame so that m v v = p ■ q. 

There is an optical theorem that relates the imaginary part of the forward Compton amplitude to 
the cross section for inelastic scattering of off-shell photons from protons. The relations precisely are 



ImtfiO,^ 2 ) = - gi{v,q 2 ) and Im H 2 {v, q 2 ) = -§■ g 2 {v, q 2 ) , (11) 



1 

v 

where g\ and g 2 are functions appearing in the cross section and are measured [18, 19, 20, 21, 22] at 
SLAC, HERMES, JLab, and elsewhere. 

Using the Compton amplitude in terms of Hi and H 2 , Eq. (10), in evaluating the inelastic part of 
the two-photon loops gives 

E 2j 



po1 ~ E F 



i 4 Q 



2am„ 

(12) 



x / {Qi + tS e Q 2 )Q 2 {(2Q 2 + Qt)Hr\*Qo,-Q 2 )-3Q 2 QlHr l Wo,-Q 2 )} , 

where we have Wick rotated the integral so that Q = —iu, Q = q, and Q 2 = Qq + Q 2 . Since are 
not measured, we obtain them from a dispersion relation, which will discussed in a subsequent section. 
Assuming no subtraction, 



7T J 2 V - V 2 



Im ^ l( ^, (13) 



where the integral is only over the inelastic region, and similarly for H 2 . 

Putting things together, neglecting m e inside the integral, and integrating what can be integrated, 
one obtains the expression 

rvrri 

A '°' = 2(TT^ (Al + A2) ' <14) 

where (with r = v 2 /Q 2 ), 

Al = \[ %{ F 2(-Q 2 )+^v [ %^)9^-Q 2 )} ^ d5) 
A 2 = -12m p f% I™ %fS 2 (j)g 2 {v,-Q 2 ). 

JO W JvtH V 

The auxiliary functions are defined by 

ft (t) = -3r + 2r 2 + 2(2 - t)^t(t + 1) = -^(t) , 



fair) = l + 2r- 2vV(t + 1). (16) 

The integral for Ai is touchy. Only the second terms comes from the proceedrue just outlined; 
it was historically thought convenient to add the first term, and then subtract the same term from the 
elastic contributions, specifically from the recoil corrections. This stratagem is what allows the electron 
mass to be taken to zero in Ai. The individual terms in Ai diverge (they would not had the electron 
mass been kept), but the whole is finite because of the Gerasimov-Drell-Hearn (GDH) [23, 24] sum 
rule, 



dv 



'imp — gi (u,0) = -K 2 , (17) 

J"th V 
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coupled with the observation that the auxiliary function /3(r) becomes unity as we approach the real 
photon point. 

The polarizability expressions have some history. A short version is that considerations of A po i 
were begun by Iddings in 1965 [12], improved by Drell and Sullivan in 1966 [13], and given in present 
notation by de Rafael in 1971 [14]. But no sufficient spin-dependent data existed, so it was several 
decades before the formula could be evaluated to a result incompatible with zero. In 2002, Faustov and 
Martynenko became the first to use gi y 2 data to obtain results inconsistent with zero [16]. They got 

A pol = (1.4 ± 0.6) ppm (18) 

However, they only used SLAC data and Ai and A 2 are sensitive to the behavior of the structure 
functions at low Q 2 . Also in 2002 there appeared analytic expressions for 512 fit to data by Simula, 
Osipenko, Ricco, and Taiuti [25], which included JLab as well as SLAC data. They did not at that 
time integrate their results to obtain A po i. Had they done so, they would have obtained A po i = (0.4 ± 
0.6) ppm [17]. 

We now have enough information to discover a bit of trouble. Table 1 summarizes how things stood 
before the 2005/2006 re-evalations of A po i . The sum of all corrections is 1 . 59 ± . 77 ppm short of what 
would be desired by experimental data. Using the Simula et al. value for A po i would make the deficit 
greater. Using other proton form factor fits (limitng ourselves to modern ones that fit the data well) in 
evaluating Az can reduce the deficit somewhat, but not by enough to ameliorate the problem [17]. 



Table 1. Corrections to hydrogenic hyperfine structure, as they could have been given in 2004. The first line with 
numbers gives the "target value" based on the experimental data and the best evaluation of the Fermi energy (8 
figures) based on known physical constants. The corrections are listed next. (The Zemach term includes a 1.53% 
correction from higher order electronic contributions [26], as well as a +0.07 ppm correction from muonic vacuum 
polarization and a +0.01 ppm correction from hadronic vacuum polarization [2].) The total of all corrections is 
1.59 ± 0.77 ppm short of the experimental value. 



Quantity 


value (ppm) 


uncertainty (ppm) 


(E hfs (e-p)/E p F ) - 1 


1 103.49 


0.01 


Aqed 


1 136.19 


0.00 


A z (using Friar & Sick [27]) 


-41.59 


0.46 




5.84 


0.15 


Apoi (from Faustov & Martyenko [16]) 


1.40 


0.60 


A p v 

weak 


0.06 




Total 


1101.90 


0.77 


Deficit 


1.59 


0.77 



The discrepancy is not large, measured in standard deviations. On the other hand, the problem is 
clearly not in statistical fluctuations of the hfs measurement one is trying to explain, so one would like to 
do better. As listed in the Table, the largest uncertainty in the corrections comes from A po i. Further, the 
polarizability corrections require knowledge of g\ and g 2 at relatively low Q 2 , and good data pressing 
farther into the required kinematic regime has relatively recently become available from JLab (the 
Thomas Jefferson National Accelerator Laboratory, in Newport News, Virginia, USA). Accordingly, 
we shall present a state-of-the-art evaluation of the polarizability correction for electronic hydrogen. To 
give away our results [17] at the outset, we essentially confirm (remarkably, given the improvements 
in data) the 2002 results of Faustov and Martynenko. 
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3. Re-evaluation of A po i 

Data for g\ (V, q 2 ) has improved due to the EG1 experiment at JLab, which ran in 2005. Data based 
on preliminary analysis is available [22]; final data is anticipated in late 2006. A sample of the new 
data is shown in Fig. 3. Since a function of two variables can be complicated to show, what is shown is 
the integral 

/•OO 

I{Q 2 ) = Am p (dv/v 2 )g x {v,-Q 2 ), (19) 

Jv th 

which differs from an integral appearing in Ai in lacking the auxiliary function. The integration was 
done by the experimenters themselves. We remind the reader that this integral is expected to reach 
the Gerasimov-Drell-Hearn value — n 2 at the real photon point, and that because of cancellations the 
difference of this integral from — k 2 is more relevant to the final answer for A po i than its absolute 
value. 




0.05 0.1 0.15 0.2 0.25 0.3 



Q 2 (GeV 2 ) 

Fig. 3. Data for (an integral over) the spin-dependent struture function gi (u, q 2 ). The data is from the EG1 exper- 
iment at Jefferson Lab from the year 2005. The Simula et al. curve is from a fit published in 2002. 



All the data shown is new; there was no polarized electron-proton scattering data available below 
Q 2 = 0.3 GeV 2 when Simula et al. and Faustov and Martynenko did their fits. A curve obtained by 
integrating Simula et a/.'s fit for g\ is also shown on the Figure; we do not have enough information to 
produce a similar curve for Faustov and Martynenko's fit. 

Integration in the region Q 2 > 0.05 GeV 2 is done using analytic fits to actual data for g\. For Q 2 
below 0.05 GeV 2 , where there is no data, we do an interpolation based on a low Q 2 expansion within 
the integral to get (with Q\ -> 0.05 GeV 2 ), 



Ai[0,Qf] = - 



dQ 2 
Q 2 



F 2 {-Q 2 )+Am % 



2 2 . 

r P K p -+- 



18m p ci 



5m 2 
4a 



■7o 



Q\ + o(Q\) . 



-Q 2 ) 



(20) 



Here rp is from the expansion of the Pauli form factor F2(— Q 2 ) = k 2 (l — rpQ 2 /6 + ...), and the 



forward spin polarizability" 

_ 2a [°° dv 
m p J Uth is 4 



7o has been evaluated from data for other purposes [28], 
[-1.01 ± 0.08 ± 0.10] x 10~ 4 fm 4 



(21) 
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The parameter c\ is defined from the slope at low Q 2 of the integral shown in Fig. 3, 



I(Q 2 ) = Am p — 9i(y, -Q 2 ) = ~4 + 8m 2 p c x Q 2 + G(Q 4 ) ; 



Jv th 



(22) 



we find and use c x = 2.95 ± 0.11 GeV 4 [17]. 

We need to comment that for A2, we need 32, and there is almost no data for 52 on the proton. 
One estimates 52 by relating it to gi using the Wandzura-Wilczek relation [29], which we shall not 
detail here. Fortunately, the auxiliary function fair) is small over the region where we need to do the 
integrals, so that even when we assigned 100% error bars to the contribution from g 2 , the effect on the 
final answer was not great. 

Our overall result is [17] 



which is similar to the Faustov-Martynenko result. This result means that the polarizability corrections 
no longer give the largest uncertainty in Table 1 . It also means that the theory deficit outlined in Table 1 
still remains, even becoming modestly larger with a smaller uncertainty limit, at (1.69 ± 0.57) ppm. 

4. Comments on the derivations of the formulas 

The polarizability corrections depend on theoretical results that are obtained using unsubtracted 
dispersion relations. One would like to know just what this means, and since there is at least a small 
discrepancy between calculation and data, one would like to be able to assess the validity of such 
dispersion relations. 

Also, the hyperfine splitting in muonic hydrogen may be measured soon. The polarizability correc- 
tions have been calculated for this case also [30], albeit only with older fits to the structure function 
data and the relevant formulas, with non-zero lepton mass everywhere, are available [30] from a single 
source, so one would like to verify these formulas. It turns out that keeping the lepton mass does not 
greatly increase calculational effort or the length of the formulas, so we can do the groundwork for the 
muonic hfs case simultaneously with the assessment of the ordinary hydrogen hfs calculation, although 
we shall not here display the formulas for non-zero lepton mass. 

The calculation begins by writing out the loop calculation using the known electron vertices and 
the definition of the Compton scattering amplitudes involving Hi and H 2 as given in Eq. (10). One 
can and should use this formalism for all the hadronic intermediate states, including the single proton 
intermediate states. The single proton intermediate states give contributions to H\ and H2 that can be 
(more-or-less) easily calculated given a photon-proton-proton vertex such as Eq. (5). For reference, we 
give the result for Hi, 



If one is calculating only the "elastic terms," one can use the result for Hf directly in the loop 
calculation and obtain a result [9] that one may call A e( . All terms in Hf will contribute, specifically 
including the F 2 term clearly visible above. A criticism of this procedure is that the proton vertex used 
is not demonstrably valid when the intermediate proton is off shell, so the above expression may or 
may not be correct overall. However, it is correct at the proton pole. 

Alternatively, one may do a unified calculation of the elastic and inelastic contributions. Since 
we don't have a direct calculation of the Hi for the inelastic case, we obtain them using dispersion 
relations. Including the elastic terms in the dispersion relation is no problem [12, 13]. One just needs 
the imaginary parts of Hf ; these are easy to obtain, and contain Dirac delta-functions that ensure the 
elastic scattering condition v = ±Q 2 / (2m p ). 



A pol = 1.3 ±0.3 ppm 



(23) 




(24) 
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Dispersion relations involve imagining one of the real variables to be a complex one an then using 
the Cauchy integral formula to find the functions Hi at a particular point in terms of an integral around 
the boundary of some region. In the present case we "disperse" in v, treating q 2 as a constant while we 
do so. Three things are needed to make the dispersion calculation work: 

• The Cauchy formula and knowing the analytic structure of the desired amplitudes. 

• The optical theorem, to relate the forward Compton Im Hi to inelastic scattering cross sections. 

• Legitimately discarding contributions from some oo contour, if the dispersion relation is to be 
"unsubtracted." 

The first two are not in question. 

For the present case, the contour of integration is illustrated in Fig. 4, where one should imagine 
the outside circle having infinite radius. The result for Hi begins its existence as 

RcsW^l , l^Wl^l, BpQs. (25) 

Ki -v 2 7T J cut v - v 2 2m J yl=00 v - V 2 



Im v 2 







y \ 













Re v 2 



Fig. 4. Contour in complex v 2 plane for applying Cauchy identity to H\ or H2. 



The numerator of the first term is the residue (Res) from the poles in v for the elastic part of Hi, 
as from Eq. (24). Note that the i 7 ^ 2 term m Hi > Eq. (24), is constant in v, certainly lacking a pole in v. 
Hence this term never enters the dispersion relation, and no F 2 term arises from the elastic contribution, 
as calculated this way. 

The second term leads to the gi term in the quantity Ai given earlier, after using the optical theorem 
to relate Im H t to gi. 

The third term is the integral over the part of the contour which is the infinite radius circle. The 
commonly quoted results for A po i, which appear in this talk, depend on dropping this term. The term 
is zero, if Hi falls to zero at infinite \v\. Assuming this is true, however, appears to be a dramatic 
assumption. It fails for Hf alone. Hence, for the assumption to succeed requires an exact cancelation 
between elastic and inelastic contributions (or a failure of Eq. (24) on the big contour). On the positive 
side are several considerations. One is that nearly the same derivation gives the GDH sum rule, which 
is checked experimentally and works, within current experimental uncertainty (8%) [31]. Also, the 
GDH sum rule has been checked in lowest order and next-to-lowest order perturbation theory in QED, 
where it appears to work [32, 33]. Finally, Regge theory suggests the Compton amplitude does fall to 
zero with energy [34], as one would like, although Regge theory famously gave wrong high v behavior 
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for spin-independent analogs of gi and g 2 [35]. Hence there are indications, though not decisive proof, 
supporting the unsubtracted dispersion relation. 

The derivation finishes by subtracting a term involving F| from the relativistic recoil term, so as 
to obtain exactly the elastic corrections A e ; = Az + A P R that were obtained (say) by Bodwin and 
Yennie for a calculation of the elastic terms only, using Eq. (5) at the photon-proton vertices and no 
dispersion theory [9]. After adding the same term to the polarizability corrections in Ai, one obtains 
the commonly quoted result for Ai [13, 14, 16]. As noted earlier, this reorganization also allows Ai 
to be finite in the m e — > limit. Beyond the historical connection, if one is comfortable with the 
unsubtracted dispersion relation, the use of the dispersion theory gives a more secure result because it 
uses only the pole part of the photon-proton-proton vertex, so that the combined elastic and inelastic 
result does not depend on the general validity of whatever photon-proton-proton vertex one uses. 

5. Conclusion 

The evaluation of the polarizability contributions to hydrogen hyperfine structure, A po i, based on 
latest proton structure function data is A po i = 1.3 ± 0.3 ppm. This is quite similar to the Faustov- 
Martyenko result, which we think is remarkable given the improvement in the data upon which it is 
based. Most of the calculated A po i comes from integration regions where the photon four-momentum 
squared is small, Q 2 < 1 GeV 2 . 

There is still a modest discrepancy between the hydrogen hfs calculation and experiment, on the 
order of 2 ppm. Where may the problem lie? It could be in the use of the unsubtracted dispersion 
relation; or it could be in the value of the Zemach radius, which taken at face value now contributes the 
largest single uncertainty among the hadronic corrections to hfs; or perhaps it is a low Q 2 surprise in 
gi or g 2 - We would like to know. It is at any rate not a statistical fluctuation in the hfs data itself. 

An interplay between the fields of atomic and nuclear or particle physics may be relevant to sorting 
out problem. For one example, the best values of the proton charge radius currently come from small 
corrections accurately measured in atomic Lamb shift [36] . The precision of the atomic measurement of 
the proton charge radius can increase markedly if the Lamb shift is measured in muonic hydrogen [37], 
which could happen in 2007, if the Paul Scherrer Institute accelerator schedule holds. In the present 
context, the charge radius is noticed by its effect on determinations of the Zemach radius. 

For ourselves, a clear continuation of the present program is to finish evaluations of the muonic 
hydrogen ground state hfs. We have the formulas with all lepton masses in place, and are currently 
waiting until the final EG1 data is released, which we think will be rather soon, before publishing with 
a numerical evaluation. 

My contributions to this subject have all been made in collaboration with Vahagn Nazaryan and 
Keith Griffioen. I thank them for the great pleasure I have had working alongside them. In addition, 
we thank Jose Goity, Savely Karshenboim, Ingo Sick, Silvano Simula, and Marc Vanderhaeghen for 
helpful discussions and information. This work was supported by the National Science Foundation 
under grants PHY-0245056 and PHY-0555600 (C.E.C.); PHY-0400332 (V.N.); and by the Department 
of Energy under contract DE-FG02-96ER41003 (K.A.G.). 

References 

1. S. G. Karshenboim, Can. J. Phys. 77, 241 (1999). 

2. A. V. Volotka, V. M. Shabaev, G. Plunien and G. Soff, Eur. Phys. J. D 33, 23 (2005). 

3. A. Dupays, A. Beswick, B. Lepetit, C. Rizzo, and D. Bakalov, Phys. Rev. A 68, 052503 (2003). 

4. S. J. Brodsky, C. E. Carlson, J. R. Hiller and D. S. Hwang, Phys. Rev. Lett. 94, 022001 (2005); Phys. 
Rev. Lett. 94, 169902 (E) (2005) [arXiv:hep-ph/0408131]. See also [5]. 

5. J. L. Friar and I. Sick, Phys. Rev. Lett. 95, 049101 (2005) [arXiv:nucl-th/0503020] and S. J. Brodsky, 
C. E. Carlson, J. R. Hiller and D. S. Hwang, Phys. Rev. Lett. 95, 049102 (2005). 

6. M. I. Eides, Phys. Rev. A 53, 2953 (1996). 

©2008 NRC Canada 



10 



unknown Vol. 99, 2008 



7. P. G. Blunden, W. Melnitchouk and J. A. Tjon, Phys. Rev. Lett. 91, 142304 (2003) [arXiv:nucl- 
th/0306076]; Y. C. Chen, A. Afanasev, S. J. Brodsky, C. E. Carlson and M. Vanderhaeghen, Phys. 
Rev. Lett. 93, 122301 (2004) [arXiv:hep-ph/0403058]; A. V. Afanasev, S. J. Brodsky, C. E. Carlson, 
Y. C. Chen and M. Vanderhaeghen, Phys. Rev. D 72, 013008 (2005) [arXiv:hep-ph/0502013]; 

J. Arrington, Phys. Rev. C 71, 015202 (2005) [arXiv:hep-ph/0408261]. 

8. C. K. Iddings and P. M. Platzman, Phys. Rev. 113, 192 (1959). 

9. G. T. Bodwin and D. R. Yennie, Phys. Rev. D 37, 498 (1988). 

10. A. P. Martynenko, Phys. Rev. A 71, 022506 (2005) [arXiv:hep-ph/0409107]. 

11. A. C. Zemach, Phys. Rev. 104, 1771 (1956). 

12. C. K. Iddings, Phys. Rev. 138, B446 (1965). 

13. S. D. Drell and J. D. Sullivan, Phys. Rev. 154, 1477 (1967). 

14. E. De Rafael, Phys. Lett. B 37, 201 (1971). 

15. P. Gnadig and J. Kuti, Phys. Lett. B 42, 241 (1972). 

16. R. N. Faustov and A. P. Martynenko, Eur. Phys. J. C 24, 281 (2002); R. N. Faustov and 
A. P. Martynenko, Phys. Atom. Nucl. 65, 265 (2002) [Yad. Fiz. 65, 291 (2002)]. 

17. V. Nazaryan, C. E. Carlson and K. A. Griffioen, Phys. Rev. Lett. 96, 163001 (2006) [arXiv:hep- 
ph/0512108]. 

18. P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 493, 19 (2000) [arXiv:hep-ph/0007248]. 

19. P. L. Anthony et al. [E155 Collaboration], Phys. Lett. B 553, 18 (2003) [arXiv:hep-ex/0204028]. 

20. R. Fatemi et al. [CLAS Collaboration], Phys. Rev. Lett. 91, 222002 (2003) [arXiv:nucl-ex/0306019], 

21. J. Yun et al, Phys. Rev. C 67 055204, (2003); S. Kuhn, private communication. 

22. A. Deur, arXiv:nucl-ex/0507022. 

23. S. B. Gerasimov, Sov. J. Nucl. Phys. 2, 430 (1966) [Yad. Fiz. 2, 598 (1966)]. 

24. S. D. Drell and A. C. Hearn, Phys. Rev. Lett. 16, 908 (1966). 

25. S. Simula, M. Osipenko, G. Ricco and M. Taiuti, Phys. Rev. D 65, 034017 (2002) [arXiv:hep- 
ph/0107036]. Silvano Simula provided us with an updated version of the code, including error estimates. 

26. S. G. Karshenboim, Phys. Lett. 225A, 97 (1997). 

27. J. L. Friar and I. Sick, Phys. Lett. B 579, 285 (2004). 

28. D. Drechsel, B. Pasquini and M. Vanderhaeghen, Phys. Rept. 378, 99 (2003) [arXiv:hep-ph/0212124]. 

29. S. Wandzura and F. Wilczek, Phys. Lett. B 72, 195 (1977). 

30. R. N. Faustov, E. V. Cherednikova and A. P. Martynenko, Nucl. Phys. A 703, 365 (2002) [arXiv:hep- 
ph/0 108044]. 

31. P. Pedroni [GDH and A2 Collaborations], AIP Conf. Proc. 814, 374 (2006); H. Dutz et al. [GDH 
Collaboration], Phys. Rev. Lett. 94, 162001 (2005). 

32. D. A. Dicus and R. Vega, Phys. Lett. B 501, 44 (2001) [arXiv:hep-ph/0011212]; 

33. G. Altarelli, N. Cabibbo and L. Maiani, Phys. Lett. B 40, 415 (1972). 

34. H. D. I. Abarbanel and S. Nussinov, Phys. Rev. 158, 1462 (1967). 

35. M. Damashek and F. J. Gilman, Phys. Rev. D 1, 1319 (1970); C. A. Dominguez, C. Ferro Fontan and 
R. Suaya, Phys. Lett. B 31, 365 (1970). 

36. P. J. Mohr and B. N. Taylor, Rev. Mod. Phys. 72, 351 (2000); and Rev. Mod. Phys. 77, 1 (2005) [2002 
COD ATA]; Cf., the electron scattering value of I. Sick, Phys. Lett. B 576, 62 (2003) [arXivmucl- 
ex/03 10008] or the spread from Sick's value to that of J. J. Kelly, Phys. Rev. C 70, 068202 (2004). 

37. A. Antognini et al, AIP Conf. Proc. 796, 253 (2005). 



©2008 NRC Canada 



